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Satellites  1961  °1,  1961  ortjl,  and  1962  Pul  are  instrumented  to 

supply  time  data  as  well  as  doppler  data.  Satellite  time  markers  are 

/  > 

generated  in  these  satellites  by  counting  down  the  satellite’s  stable 
oscillator  and  producing  a  unique  modulating  on  the  doppler  carriers  each 
time  a  fixed  number  of  oscillator  cycles  have  occurred.  The  TRANET  doppler 
receiving  stations  are  equipped  to  recognize  this  modulation  and  to 
automatically  punch  (on  the  same  paper  tape  that  doppler  data  is  punched) 
the  station's  time  (WWV)  at  which  this  special  modulation  is  received. 
After  appropriate  calibration,  this  time  data  can  serve  as  a  secondary 
time  standard.  This  report  presents  the  information  required  to  recognize 
and  process  satellite  time  data  and  to  determine  the  WWV  time  at  which 
the  satellites  transmitted  satellite  time  markers. 

The;  report  is  divided  into  two  sections,  first  section  presents 


The;  report  is  divided  into  two  sections.  first  section  presents 

the  basic  equations  that  connect  the  satellite  transmitter  frequency,  the 
station’s  satellite  time  data,  and  the  times  that  the  satellites  transmitted 
the  time  markers.  The  second  section  presents  the  information  required  to 
recognize,  re-format,  and  reduce  the  experimental  satellite  time  data  to  a 
form  appropriate  to  the  theory  presented  in  the  first  section. 


£ 

In  this  report,  WWV  Time  is  used  to  designate  WWV’ s  estima 
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I.  DERIVATION  OP  EQUATIONS  FOR  SATELLITE  TIME 

Tills  section  presents  the  equations  connecting  the  satellite's 
oscillator  frequency,  the  stations'  data  when  the  satellite  transmitted 
time  markers,  and  inferred  values  for  the  times  when  time  markers  are 
transmitted  relative  to  a  chosen  standard  of  time.  The  section  is 
divided  into  three  sub-sections.  The  first  sub-section  presents  the 
notation  for  important  quantities  which  are  used  throughout  this  report. 
The  second  sub- section  presents  the  equations  to  determine  the  absolute 
frequency  of  the  satellite  as  a  function  of  time.  The  third  sub-section 
presents  one  method  for  statistically  inferring  when  the  satellite 
transmitted  time  markers.  Also  included  in  the  third  section  is  a 
suggestion  for  a  simplified  equation  to  predict  into  the  future  the  times 
that  satellite  time  markers  are  transmitted  for  station  alert  purposes. 

A.  Notation 

* 

1.  Satellite  Frequency 

fR  =  satellite  reference  frequency  (mHz); 

=  standard  offset  of  satellite  frequency  (hz/mHz); 

fg  0  =  standard  satellite  transmitter  frequency  (hz), 

=  fR(l06  +  Aus); 


1  az  -  1  cps . ,  and  1  mHz  =  1  me . 
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f„(t)  =  estimated  satellite  freqxiency  as  function  of 
°  titue  (Hz); 


5fg(t)  =  estimated  deviation  of  satellite  frequency 

from  standard  (Hz), 


ig(t)  -  ig^Q; 


Ctig(t »t^ ) 3  =  functions  chosen  for  expansion  of  Bfg(t), 


Bfs(t)  =  E  a  £  «A(t,t  ); 

1=1  1 

=  epoch  of  satellite  frequency  equation. 


(a„  ^ }  =  expansion  coefficients  of  deviation  of  satellite 
i''  frequency  for  the  {u^(t,tQ)}j 

T> 

vl(t^t0)  =  l  dt’ 

X0 

5f  =  experimental  determination  of  the  deviation  of  the 

satellite  frequency  from  its  standard,  scaled  to  fR, 
for  the  y-th  pass, 


cy 


time  of  closest  approach  of  the  satellite  to  receiving 
station  for  the  Y-th  pass,  (sec). 


2.  Satellite  Time 


{  £  1 

tg  =  epoch  of  satellite  time  equation; 

ZTg(t)  =  number  of  cycles  of  satellite  oscillator  (scaled  to  f^) 
that  have  occurred  after  tg^5'; 

g  =  number  of  cycles  of  satellite  oscillator  (scaled  to  f^) 
'  that  occur  between  two  successive  transmissions  of 
time  markers; 
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K  =  number  of  satellite  time  markers  that  have  heen 


transmitted  since  \ s  ^ ; 


"fell 

tj.  =  inferi*ed  time  of  K  time  marker  transmission; 
fs )  th 

t.  '  ^  =  experimentally  measured  time  of  K  time  marker 
transmission; 


Int.  [X(t)l  =  rounded  integer  part  of  the  nimiher  X(t). 


B.  Inferred  Satellite  Frequency  as  a  Function  of  Time 

Within  a  few  days  after  successful  launch  of  a  satellite  containing 
stable  oscillators,  a  standard  offset  of  the  oscillator  frequency  from 
exactly  one  h5s  is  chosen.  This  standard  offset,  Aug,  normally  remains 
fixed  for  the  useful  life  of  the  satellite.  Table  I  gives  the  value  of 
ius  for  those  two  satellites  currently  in  orbit  which  transmit  time 
markers . 


Satellite  (hz/mHz 

19bl  ol  -35-333333 

1962  Sul  -77.000000 

TABLE  I.  Standard  Offset  of  Satellite  Oscillators 

An  accurate  experimental  measurement  of  the  satellite  oscillator  frequency 
can  be  obtained  during  pre-processing  of  the  doppler  data  by  adjusting  the 
receiving  station's  position  to  provide  the  best  fit  to  the  experimental 
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(l) 

doppler  data . v  *  This  experimental  measurement  of  the  satellite  oscillator 
frequency  for  the  Y-th  pass  is  usually  scaled  to  some  convenient  reference 
frequency,  f^,  and  is  given  as  a  correction  term,  &f  ,  to  the  standard 

/f 

satellite  freqxteney,  fg  Q  =  f^(lO°  +  Aug). 

Consequently,  letting  the  time  of  closest  approach  for  the 

Y-th  nsss.  t  ,  be  the  time  associated  with  the  experimental  measurement 
CY 

of  frequency,  &f  ,  it  can  be  seen  that  the  deviation  of  the  satellite 

frequency  from  its  standard  frequency,  &fg(t),  can  be  directly  inferred 

from  the  experimental  data,  C&f^} .  Let  the  deviation  of  the  satellite 

frequency  be  parameterized  in  the  form: 

n_ 

1 

&fs(t)  =  £af^  u£(t,t0)  (1.1) 

4=1 

where  the  u^(t,tQ)  are  known  functions  of  the  time  and  chosen  to 
accurately  represent  the  variations  of  the  satellite  oscillator 

* 

frequency  as  a  function  of  the  time  t,  and  the  frequency  epoch,  t^. 

Choosing  a  convenient  epoch;  t^,  for  equation  (l.l),  a  least  squares 
determination  for  the  expansion  coefficients,  [a^ 


*  1961  oTJl  data  indicates  that  for  active  temperature  controlled  satellite 

oscillators,  ir,  =  1  and  (t,t^)  =  t  -  are  adequate  for  parameterizing 
the  variations  of  the  satellite  frequency,  so  long  as  the  ambient  satellite 
temperature  remains  within  the  control  range. 

(l) "Orbit  Improvement  Program",  APL/JHU  Report  TG-401. 
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F(af)  =  n 


n 

1  J? 


S  w. 


P  Y=1 


4  *t,i  VVV  '  8fV 


(1.2) 


2F(a  ) 
5a.  « 


=  0;  l=l,2,....n 


s  •  •  •  • 


where  n 

P 


total  number  of  satellite  passes, 


provides  a  method,  for  inferring  the  satellite  oscillator  frequency  at 
any  desired  time.  “Hie  statistical  weighting  factors,  w^,  should  be 
chosen,  to  correspond  to  an  a  priori  estimate  of  the  quality  of  the 
doppler  data  ana  the  accuracy  of  the  satellite  orbit  used  to  infer  the 
experimental  {§f^} . 


C.  Inferred  Times  when  the  Satellite  Transmits  Time  Markers 

The  satellites  considered  in  this  report  initiate  the  transmissions 
of  time  markers  by  counting  the  number  of  cycles  generated  by  the  satellite 
oscillator  (positive  zero  crossings  of  the  satellite  oscillator  voltage) 
and  triggering  the  transmission  each  time  that  the  count  reaches  a 
predetermined  integer  number.  Consequently,  except  for  a  scaling  factor, 
the  oscillator  frequency  represents  the  satellite  'clock  rate'  and  all 
that  is  required  to  determine  accurately  the  times  that  markers  are 
transmitted  is  to,  l)  determine  accurately  the  satellite  frequency,  and 
2)  determine  one  time  of  transmission  of  a  marker  relative  to  some  accepted 
standard  of  time  (the  epoch  of  the  time  equation).  This  section  presents 
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the  information  required  to  determine  the  scaling  factor  for  the  clock 
rate  and  the  equations  to  determine  the  time  of  transmission  of  one 
time  marker,  given  experimental  measurements  of  the  times  of  transmission 
of  time  markers  relative  to  an  accepted  time  standard.  In  addition, 
this  section  pie  sent s  the  equations  necessary  to  check  the  station 
clock  of  any  receiving  station.  Finally,  this  section  presents  one 
fora  of  a  simplified  equation  that  is  convenient  to  use  for  extrapolating 
marker  transmission  times  into  the  future  to  an  accuracy  sufficient  for 
station  alert  purposes. 

If  3f_(t)  is  the  number  of  cycles  that  the  satellite  oscillator 

O 

(scaled  to  f.,)  has  generated  since  some  epoch,  then  the  change  in  Ng(t) 
over  a  short  time  interval  can  be  used  as  a  measure  of  the  satellite 
frequency  with  sufficient  accuracy. 

ft  (t  +  lit)  -  N  (t) 

— - it - 5 - fs<*)  -  fs,o  +  Ms(t),  d.3) 


or,  zo  sufficient  accuracy,  NQ(t)  can  be  represented  by: 


HS <*>  -  *  fs,o(t  -  *o>  *  1  at’  Ms(t'). 


(1.1.) 


let  E_  c  be  the  number  of  cycles  generated  by  the  oscillator  between 
two  successive  time  marker  transmissions,  (again  scaled  to  f^)  and  let 

t 


Hs(*)  n  i 

*B(t)  “  =  W  + 


j dt*  8fs(t*) 


(1.5) 
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Finally,  if  the  epoch  for  Kg(t)  is  the  time,  tQS\  that  a  marker  was 
transmitted,  the  time  that  the  Kth  time  marker  is  transmitted,  tj ,  is 
determined  by  tlie  equation: 


Integer  K  =  Kg(t^) 


!S(t0} 


s,o 


N. 


C,S 


(tk  '  t0)  +  N 


c,s 


dt'  5f_(t' ),  (l.6A) 


with  the  boundary  condition  that: 


0  =  Kg(t0) 


S  .0  fj_{s) 


W 


C,S 


(C'  -  + 


0  o'  w, 


c,s 


■  (s) 


dt'  Sfg(t*). 


(1.6b) 


Combining  equations  (1.6), 


rr  S.O  ^(S)x 

K  "ITT  --  '  °o  } 

v  .O 


H, 


C.S 


“k 


.(b) 


j  dt '  6fs(t') 
feo 


dt»  5fg(t') 


(1-7) 


where  K  is  always  an  integer. 


Finally,  substituting  equation  (l.l)  of  the  previous  sub-section 
into  equation  (l.7),and  letting: 


v^(t,t0)  =  j  dt'  u^(t,t0). 


(1.8A) 
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then 


S.Q 


M, 


C,S 


(\  -  4S)> * 


N, 


c,s 


f 

j£x  af,i[v^Vto)  -  ^o^V1 


(1 


Table  I  in  Section  I. A.  listed  \?alues  of  AUg  from  which  fg  Q  can  he 

calculated.  Table  II.  giran  he low ,  lists  the  values  of  Ng  g  for  the  three 

satellites  of  interest  in  this  report.  Assuming  that  the  frequency  epoch, 

t  ,  has  been  chosen  and  the  values  of  the  frequency  expansion  coefficients, 

u 

a_  ,, ,  have  been  inferred  from  doppler  data  as  discussed  in  section  I.  A,, 

M  (s 

there  remain  in  equations  (1.8)  only  three  possible  unknowns :  K,  t^,  or  t^ 

It  can  be  seen  that  if  any  two  of  these  three  are  known,  the  third  can  be 

found  from  the  above  equations.  Three  cases  are  now  discussed,  the  case 

number  depending  upon  which  of  the  three  variables  is  considered  unknown. 


Sav,e_Lli  L»e 
1961  ol 

1961  cr  1 

1962  Sul 


Nc  gat  3  mHz. 

100712 448 

67092480 

268435456 


TABiE  11.  Number  of  Oscillator  Cycles  Occurring  between  Time  Markers 


Case  1:  The  Number  of  the  Time  Marker,  K,  is  Unknown. 

From  the  numbers  given  in  Tables  I  and  II,  it  can  be  seen  that 
two  time  markers  can  not  occur  more  frequently  than  about  every  20  seconds. 
Consequently,  given  an  approximate  time  of  transmission  of  a  time  marker, 

(s) 

ty,  and  an  approximate  value  for  the  time  of  the  zeroth  time  marker 


TM  Wvm  HopAm\  Uoivenity 

APPUt*  PMYAtCt  LAMA  ATOPY 

Spriog.,  Wanf-l*** 

-  10  - 


transmission,  (epoch  for  time  equation),  both  accurate  to  within 

(e) 

a  few  seconds,  the  mine  of  the  integer  K  corresponding  to  t>  is 
uniquely  given  by: 


K  =  lot 


W 


j0_/t(E)  _  t(s)}  +  .  1  yf  a  (v  (t(E)  t  1  v  (t^  t  )) 

„  ''  k  0  '  N„  „  A  VlAAk  >V  Yl[to  o )} 


C,S 


C,S  1=1 


(1.9) 


where  Int  [  J  represents  the  rounded  integer  part  of  the  quantity  contained 
in  the  bracket. 


Case  2:  "Hie  inferred  Time,  t^,  is 


Unknown . 


.(s) 


Given  an  accurate  value  of  the  epoch,  t^  ,  and  an  approximate 

(E) 

value  for  the  time  of  transmission  of  a  time  marker,  t^  ,  the  corresponding 

integer  S  can  be  found  from  equation  (l-9)-  Then,  the  inferred  value,  t^, 

(E) 

which  lies  closest  to  t,_  '  can  be  found  from  the  equation: 


,  _,(s) 


N 


C-  =  U 


.  K  - 

s.o 


7“"  zl±  af,i  ^(tk,tQ)  ’  vi(toS)^o)]  (1 


This  is  a  transcendental  equation  in  the  unknown  t,  .  However,  the 

k  ’ 

contribution  of  the  sum  to  equation  (l.lOA)  is  so  small  that  an  iterative 

(m) 

procedure  converges  very  rapidly  (about  two  iterations).  Let  '  be 
the  result  after  m  iterations.  Then: 


n. 


Jm+1)  As)  hC,S  .  _ 1_  f 

f  ^  f 

~S,0  S,0  1=1 


.10, 


(1.10B) 
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where 


.(0)  _  t(E) 

'k  ~  \  * 


(s) 

Case  3:  Hie  Epoch,  ' ,  is  Unknown 

f  s  } 

The  solution  for  t^  is  the  most  laborious  of  the  three  possible 

cases.  Since,  given  experimental  data  on  the  times  of  marker  transmissions , 

(s) 

it  is  most  important  to  obtain  an  accurate  estimate  of  the  epoch,  t^  , 
this  case  is  discussed  in  more  detail.  The  detailed  procedures  presented 
here  are  those  that  have  been  programmed  by  The  Applied  Physics  Laboratory 
and  have  proved  to  be  adequate. 


Clearly,  if  an  approximate  value  for  the  epoch  is  known,  an 


iterative  procedure  analogous  to  the  one  given  for  Case  2  can  be  used. 

(s) 

furthermore ,  even  if  the  starting  value  for  '  in  the  iterative  procedure 

is  grossly  in  error,  the  contribution  of  the  sums  over  the  functions, 

{v,(t,t  )},  is  so  small  that  the  iterative  procedure  still  converges 
x>  'O' 

rapidly.  Because  of  this,  it  has  been  found  that  obtaining  an  accurate 

(s) 

starting  value  for  ,  so  that  the  iterative  procedure  is  efficient,  is 
not  she  principal  problem.  The  greatest  difficulty  is  in  determining 


which  data  is  spurious  so  that  only  valid  data  is  used  in  the  final 
cosspurauicn  for  t Because  of  this,  the  obvious  procedure  of  choosing 
the  epoch  to  be  near  one  of  the  experimental  data  points  does  not 
necessarily  decrease  computing  time. 

(s) 

it  is  very  convenient  to  choose  the  epoch,  ,  to  be  close 

to  the  epoch  for  the  frequency  equation.  In  the  APL  programs,  the  epoch 
for  the  frequency  equation  is  chosen  to  he  identical  to  the  epoch  for 


Tfcfc  Jdhru  HepAtSm  UnWttvtV 

APPLIED  PHY*  ICt  LA  MR  ATOPY 

S!Si«r  Sptw^  W**vl*>d 

-  12  - 


the  satellite  orbit.  Using  this  epoch  as  a  starting  value  for  the 

(s\ 

procedure  to  determine  t^  1  insures  that  all  epochs  remain  near  each 
other  in  time  -  an  obvious  advantage  from  the  standpoint  of  "bookkeeping". 
In  detail,  the  APL  prescription  for  determining  the  epoch, 

(s) 

ts  ‘  is  the  following.  Given , 

o 


t(s) 
\  * 


k  =  1,  2 


=  experimental  time  data, 


w  ,  k  =  1,  2,  ....  n^,  =  a  priori  statistical  wts, 


=  epoch  for  frequency  equation. 


Set: 


(s,0)  _  j. 


0 


where  is  the  m-th  iterative  result  for  tpS .  For  every  data  point. 


CGSSXJUti© 


(0)  fs,0  rJE)  (s,0) 

^  -  Nc>s  K  -  *0 


c,s 


nf 

E 

1=1 


af,i  Cv^(4E)^ o}  -  W  V1' 


]. 


(1.11A) 
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Then,  for  every  data  point,  compute: 


W 

4S,1)  -  ‘o'*0) +  T§f-  kf,  \<4  -  KS) 


"t 


-  X  "fc^o8'1’  -  ‘o8’0*  -  <40)  -  40)>]2- 


k=l 


.(s) 


This  equation  provides  a  good  starting  value  for  t^  ,  providing  that 
a  sga 1 1  percentage  of  the  time  data  points  are  spurious.  The  normal 
iteration  loop  is  given  as  follows. 

Compute  just  once  mol's,  for  each  data  point. 


\ 


=  Int 


"s.o  /.(E)  .(s,0)s 

iNc,s  (°k  "  °° 


—  2  ^  tn)  -  v„(t^Sil),tn)) 


"^cs  f=i  f>A  z  k  5  0  1  0 


k  -  1,  2,  ...  .ru,. 


and  then,  iterate  with  each  fixed, 


,.(s,m)  _  .(E)  _  . (s ,m)  _  kC,S 

^5.  -P 

S,0 


O.k  ~k  "0  t„  „  x. 


iTT  /  -  u<tos’m'’to)>’ 

S;0 


(s,m) 


(1.11B) 

only 


(1.12A) 


(l.lSB^ 
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a 


"T 

E  w. 


■At. 


At(s,m) 

#  \  ,  ,  k  O.k 

_  k=l _ * 

'0 


“t 

e'  w. 


k=l 


(l.lSBg) 


ffirl 


“t 


k=l 


'  k 


“T 

e' 

k=l 


(1.12B  ) 


where  the  primed  sum  includes  only  those  data  points  for  which 


£  Q 


m 


(1.12C) 


and 


^  (s,s+l)  _  ^(s,m)  A+(s,m) 

°0  “  °0  zo 


(1.12D) 


Breakout  of  the  iteration  loop  occurs  when  At^S,m  becomes  less  than  10  ^ 
seconds  ana  the  number  of  valid  data  points  has  remained  constant  for  two 
iterations . 

(s) 

saving  inferred  a  value  for  the  epoch,  t^  ,  any  experimental 

data  point  can  be  checked  by  using  the  procedures  shown  in  Cases  1  and  2. 

Clearly,  only  data  originating  from  'time  standard  stations'  should  be 

(s) 

used  in  evaluating  the  epoch  t^  ' . 
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The  procedures  of  Case  2  can  be  used  to  compute  an  ephemeris  of 
times  when  time  markers  are  transmitted  and  can  be  used  to  extrapolate 
time  marker  transmission  times  into  the  future.  However,  for  transmittal 
of  time  ephemeris  data  from  one  agency  to  another,  the  following  simplified 
equation  for  providing  time  marker  transmission  times  is  convenient. 

Tiie  form  of  the  equation  is  chosen  to  be 

tk  -  t<*>  +  c,  K  +  C2  K2 

where  the  constants  C1  and  are  chosen  such  that  the  t  as  computed 
by  this  equation  agrees  to  sufficient  accuracy  with  the  more  accurate 
equations  given  previously.  Experimental  data  from  1961  oT|l  indicates 
that  for  oscillators  under  active  temperature  control,  an  equation  of 
this  form  will  remain  valid  for  over  one  week  to  an  accuracy  greater 
than  10  milliseconds.  Consequently,  such  an  equation  is  ideal  for 
station  alerting  purposes. 


II.  RECOGNITION  AND  PRE-PROCESSING  OF  SATELLITE  TIME  DATA 

Satellite  time  data  for  any  one  satellite  pass,  as  punched  on 
paper  tape  and  received  at  the  Satellite  Control  Center,  is  interspersed 
with  the  doppler  data  for  the  same  satellite  pass.  This  section  presents 
the  information  required  to  differentiate  between  doppler  and  time  data 

points,  and  to  reduce  the  experimental  time  data  to  experimental  times, 

(E)  * 

,  when  the  appropriate  satellite  transmitted  satellite  time  markers. 

This  report  does  not  consider  the  change  in  data  format  necessitated  when 
7  digits  instead  of  6  are  used  for  the  doppler  period  count. 
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A  typical  doppler  data  point,  when  pun  died  on  TWX  paper  tape, 
lias  the  format: 

CL  XXXXX  YYYYYY 

where 

C  =  TWX  carriage  return 
L  =  TWX  line- feed 
XXXXX  =  t 

u 

=  First  5  TWX  cliaracters  (integers)  denoting  the 
station  time  of  the  p-th  doppler  data  point,  (sec). 

YYYYYY  =  Tp 

=  Last  6  TWX  cliaracters  (integers)  denoting  the 
doppler  period  count,  (microseconds). 

when  a  station  is  taking  satellite  time  data,  any  arbitrary 
doppler  data  point  can  be  suppressed  and  replaced  by  a  satellite  time 
data  point.  He  format  of  the  satellite  time  data  point  can  be  in 
either  of  two  different  formats,  depending  upon  the  type  of  timing 
equipment  installed  in  the  static.  For  identification  purposes,  these 
two  formats  will  be  denoted  as  APL  or  NOTS  format.  They  are: 

APL  Format 
CL  9ZZZZ  XXXXXO 

where 

C  =  TWX  carriage  return 
L  =  TWX  line-feed 
XXXXX  =  t 

P 


NOTS  Format 
CL  XXXXX  OOZZZZ 
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( Format  cant'd.) 

=  5  TWX  characters  (integers)  giving  the  integer  part 
when  the  satellite  time  mark  is  received,  (sec). 

1MIL  =  AT 

=  4  TWX  diaracters  (integers)  giving  one  correction 
to  he  applied  to  t  (tenths  of  milliseconds). 

p, 

With  either  format  a  satellite  time  data  point  can  he  uniquely 
identified  when  interspersed  with  doppler  data  because  the  TWX  9  in  the 
APL  format  and  the  first  TWX  0  (zero)  in  the  NOTS  format  are  illegal 
characters  for  a  doppler  point.  The  following  flow  chart  is  an 
example  of  a  testing  network  to  correctly  identify  each  data  point. 

Upon  identification  of  a  satellite  time  data  point,  the 
TWX  C  and  L  should  be  deleted,  the  remaining  11  TWX  digits  appropriately 
unpacked,  and  transformed  from  BCI  to  the  two  floating  point  numbers. 

tM  =  integer  part  of  satellite  time  in  seconds 

AT  =  correction  for  satellite  time  in  units  of  10  ^ 
u 

seconds. 

The  reconstruction  of  the  time  (WWV)  that  the  satellite 
transmitted  a  marker  is  given  below.  'This  reconstruction  also  depends 
upon  the  type  of  station  equipment  (APL  or  NOTS)  and,  if  obtained  from 
APL  type  equipment,  further  depends  upon  the  satellite  transmitting 
the  time  marker. 
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c 


A.  Reconstruction  of  APL  Type  Time. 

1  ' Iv  ‘  ‘ "■  ■ ,i'’’  *  ** '  ~ 

Where 

] 

time  (station's  estimate  of  WWV  time)  at  which  time 
markei'  transmitted  (sec), 
defined  above, 

delay  time  in  station  time  gear  for  the  s-th 
satellite  (sec), 

station  clock  correction  given  in  data  header  (sec), 
transmission  time  from  satellite  to  station  (sec), 

and  where; 


.(E)  _ 


t  .At  = 

U'  p, 

,(s)  _ 

°D  - 


Atw  = 


At 


TR 


Satellite 


Value  of  tp  (sec) 


1961  ol 
1961  aT,l 


3.27692  x  10' 
1.66406  x  10 


1962  Bpl 


0.0 


B.  Reconstruction  of  N0TS  Type  Time 


4E>  *  [V  +  1  -  10'h  -  4V  +  4tn  -  “m 


(2.2) 
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where 


.(E) 


t  .  AT 
P'  u 

At_ 


ioD 


At 


TR 


=  time  (station’s  estimate  of  WWV  time)  at  which  time 
marker  transmitted  by  satellite,  (sec), 

=  defined  above, 

=  station  clock  correction  given  in  data  header,  (sec), 
=  WOTS  station  delay  time 

=  ti'ansmission  time  from  satellite  to  station,  (sec). 


(e  ) 

Equations  (2.l)  and  (2.2)  define  t^.  '  which  is  the  stations  best  estimate 
of  the  time  (WWV)  that  a  satellite  time  mark  is  transmitted  from  the 
satellite.  Typically,  from  one  to  fifty  such  times  can  be  received 
during  a  single  pass.  Currently  t^  for  the  NOTS  type  time  is  zero. 
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